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Wc report on a search for the pair production of second generation scalar leptoquarks {LQ) 
in pp collisions at the center of mass energy ^/s = 1.96 TeV using a data set corresponding to an 
integrated luminosity of 1.0 fb^^ collected with the DO experiment at the Fcrmilab Tevatron Collider. 
Topologies arising from the LQLQ —^ jiqvq and LQLQ —^ nqnq decay modes are investigated. No 
excess of data over the standard model prediction is observed and upper limits on the leptoquaxk 
pair production cross section axe derived at the 95% C.L. as a function of the leptoquark mass 
and the branching fraction (3 for the decay LQ —> jjtq. These are interpreted as lower limits on 
the leptoquark mass as a function of 3. For 13 = 1 (0.5), scalar second generation leptoquarks with 
masses up to 316 GeV (270 GeV) are excluded. 
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PACS numbers: 14.80.-j, 13.85.Rm 



The observed symmetry between lepton and quark 
generations could be explained by new gauge bosons in- 
troducing couplings between the lepton and quark sec- 
tors. Such particles, commonly referred to as lepto- 
quarks [H , would carry both lepton and baryon quantum 
numbers as well as fractional electric charge. Extensions 
of the standard model (SM) based on a larger gauge sym- 
metry group usually predict the existence of massive lep- 
toquarks. Experimental bounds on lepton number viola- 
tion, proton decay, and flavor changing neutral currents 
constrain hypothetical leptoquarks with masses of several 
hundred GeV to a few TeV to couple only to one quark 
and one lepton family, via processes conserving both lep- 
ton and baryon quantum numbers. Three generations of 
leptoquarks can thus be distinguished by considering the 
lepton observed in the leptoquark decay. 

At the Fermilab Tevatron pp Collider, leptoquarks 
would predominantly be produced in pairs via qq an- 
nihilation into a gluon in the s-channel independently 
of the unknown coupling A between the leptoquark and 
its associated lepton and quark. Thus, for scalar lep- 
toquarks the production cross section depends only on 
the strong coupling constant and the assumed leptoquark 
mass. The additional contribution to leptoquark pair 
production from t-channel lepton exchange with a cross 
section proportional to A'^ can be neglected. When adopt- 
ing the assumption that the leptoquarks couple to leptons 
and quarks of the same generation, the i-channel process 
is further suppressed for second and third generation lep- 
toquarks due to the vanishing parton distribution func- 
tions (PDFs) at high proton momentum fractions x for 
partons other than u and d quarks. Leptoquarks can 
decay either into a charged lepton and a quark with a 
branching fraction /3 or into a neutrino and quark with 
a branching fraction (1 — /3), assuming that the lepto- 
quark mass is much larger than the masses of its decay 
products, which is generally the case for first and second 
generation leptoquarks. Consequently leptoquark pair 
production could lead to three characteristic final states: 
£~^q£~q, t^qvq^ and vqvq^ with branching fractions /3^, 
2/3(1 — /3), and (1 — respectively. 

This Letter describes a search for second generation 
scalar leptoquark pair production in the decay modes 
LQLQ — !■ iiqvq and LQLQ — > fiqfiq using pp collisions 
at the center of mass energy = 1.96 TeV recorded with 
the DO detector. These channels lead to topologies with 
one muon, missing transverse energy (from which the 
transverse momentum of the neutrino is inferred), and 
two jets {fi^rjj signature), or with two muons and two 
jets {nnjj signature), respectively. At Run II of the Teva- 
tron, the dimuon signal was previously studied by the 
DO collaboration 0] while the CDF collaboration stud- 
ied both dimuon and single muon signals with smaller 



data sets. Since one of the muons of the LQLQ fJ.qfJ-q 
decay mode might not be reconstructed, this signal con- 
tributes to the single muon signature as well, which is 
taken into account in our analysis. The contribution of 
LQLQ iJ-qi^q in the dimuon selection can be neglected 
due to the small probability for a jet to mimic an isolated 
muon. 

The signal sensitivity for both signatures depends on (3. 
The branching fractions for the decay modes LQLQ — > 
fj.qfj.q and LQLQ — > fiqvq are maximal at /3 = 1 and 
P = 0.5, respectively. The branching fraction of both 
decay modes vanishes for /3 = 0, where the leptoquark 
would decay exclusively into a neutrino and quark. The 
resulting acoplanar jet topology has not been investi- 
gated in this analysis, but was studied recently using a 
larger data set [J]. 

The DO detector [1] is designed to maximize the de- 
tection and identification of particles arising from pp in- 
teractions and is constructed of dedicated subsystems ar- 
ranged around the interaction point. The central track- 
ing system, located at the innermost part of the detec- 
tor, consists of a silicon microstrip tracker (SMT) and 
a central fiber tracker (CFT) which cover the pseudora- 
pidity regions |?7| < 3 and |7/| < 2.5, respectively. The 
pseudorapidity is defined as = — ln[tan(0/2)] where 
6 is the polar angle with the proton beam. A 2 T su- 
perconducting solenoidal magnet is positioned between 
the tracking system and the surrounding central and for- 
ward preshower detectors. The calorimeter is composed 
of a central section (CC) which covers \ri\ < 1.1 and two 
end calorimeters (EC) that extend coverage to \ri\ « 4.2. 
The three calorimeter sections are housed in their own 
cryostats and consist of successive layers of mostly ura- 
nium absorbers and active liquid argon The muon 
system [Tj] is located outside the calorimeter and covers 
the region jr;! < 2. It consists of a layer of drift tubes 
and scintillation counters before 1.8 T iron toroids and 
two similar layers outside the magnets. 

This search for leptoquark pair production is based on 
an integrated luminosity of 1.0 fb~^. The data samples 
for the single muon and dimuon analyses are selected 
with combinations of single muon triggers jsl . 

Muons are identified in the region jryj < 2 using track 
segments found in the muon detector which are required 
to have hits in both the drift tubes and the scintillation 
counters. The segments are matched to tracks recon- 
structed in the central tracking system which determine 
the muon momenta. Muons with a transverse momen- 
tum pt > 20 GeV are kept. A veto on cosmic ray muons 
is applied which is based on timing information in the 
muon system and the removal of events with an appar- 
ent muon pair back-to-back in pseudorapidity. For the 
fi^xjj selection, a tight muon identification is applied in 
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order to suppress events without prompt muons. Tight 
quality is defined by requiring additional hits in the muon 
detector and that the muon track be isolated from other 
tracks with the sum of the transverse momenta of all 
other tracks in a cone defined in terms of 77 and azimuth 
(f) with radius TZ = y%Ar^p + (A0)^ < 0.5 around the 
muon less than 2.5 GeV [81]. The muon isolation is fur- 
ther improved by selecting events in which the energy 
measured in the calorimeter in a hollow cone of radius 
0.1 < TZ < 0.4 around the muon is less than 2.5 GeV. 
In the case of the /Lt/ijj selection, where the existence of 
two reconstructed muons allows for a better separation 
of signal and background, the hit requirement is loos- 
ened and only a track isolation criterion with a threshold 
of 4 GeV on the transverse momentum sum is required 
(loose quality). 

Jets are reconstructed with an iterative, midpoint cone 
algorithm with a cone radius of 0.5 Only jets found 
within I77I < 2.5 and with px > 25 GeV are kept. The jet 
energies are calibrated as a function of the jet transverse 
energy and 77 . 

The transverse momentum of a final state neutrino can 
be inferred from the -^t, calculated as the vector sum of 
the transverse energies in the calorimeter cells which is 
corrected with the transverse momenta of the selected 
muons and the jet energy calibration. 

The main SM background to the pair production of 
leptoquarks followed by their decay into the fiqvq and 
final states is the associated production of jets with 
W and Z/j* bosons, respectively. Vector boson produc- 



tion is simulated using the ALPGEN event generator [11 1 
which is interfaced to PYTHIA [i3] for the simulation of 
parton showering and hadronization. Samples with up to 
five or three partons in addition to the W or Z/j* bo- 
son, respectively, are generated and combined using the 
MLM matching prescription [l3|. For these samples as 
well as for the ones described below, the CTEQ6L1 PDF 
sets [l3l are used. Additional samples for W and Z/j* 
boson production are generated with the PYTHIA event 
generator and utilized to determine uncertainties in the 
transverse momentum shape of the associated jets. Top 
quark pair production also contributes to the analyzed 
final states and is simulated with PYTHIA assuming a 
top quark mass of 175 GeV. The multijet background, 
in which a jet is misidentified as an isolated muon, is 
estimated from data for the fi^xjj selection by invert- 
ing the muon isolation requirement and normalizing the 
obtained rate to data with standard isolation in a re- 
gion with < 10 GeV, which is dominated by multijet 
events. In the fJ.njj selection, the requirement of two re- 
constructed muons with large invariant mass suppresses 
background from multijet production. The background 
contribution from diboson or single top quark produc- 
tion, which could potentially contribute in both single 
muon and dimuon analyses, are found to be negligible 
as well. Leptoquark pair production is simulated using 



PYTHIA for leptoquark masses ranging from 140 GeV to 
320 GeV, corresponding to production cross sections be- 
tween 2.4 pb and 0.0074 pb. The generated events are 
processed through a full simulation of the DO detector 
based on GEANT [l^ ■ 

Leptoquark candidate events are selected by requir- 
ing at least two jets with pT > 25 GeV and muons with 
Pt > 20 GeV. For the nfijj signature at least two muons 
of loose quality are required, while for the ^^rjj selec- 
tion, events with exactly one muon of tight quality are 
kept and a veto on events with additional loose muons is 
applied to ensure that the samples have no overlap. 

For the /x/Ltjj sample, the dimuon invariant mass 
M{fi, /z) reconstructed from the two leading muons (i.e. 
muons with highest px) is required to be larger than 
50 GeV. The numbers of Z/j* and W boson events are 
simultaneously normalized to data with a common scale 
factor in the region of the Z boson resonance (defined 
as 60 GeV < M{^i,^j.) < 120 GeV) after aU the preced- 
ing cuts. At this stage, 913 data events and 930 ± 151 
expected background events remain, with a signal effi- 
ciency of 39.7% for an assumed leptoquark mass M^q — 
280 GeV. 

To account for the muon px resolution which can result 
in an overestimation of M(fi,^) and a mismeasurement 
of and in order to enhance the separation between 
signal and background, the reconstructed dimuon invari- 
ant mass is corrected by substracting the projection of 
onto the leading muon direction from the leading 
muon Pt- In this way is minimized, which is consis- 
tent with the expectation that there is no genuine hi 
both signal and Z/j* background. The minimum of the 
initial and corrected dimuon invariant mass M(/z, /z)™'" 
is utilized as a selection variable to improve the rejec- 
tion of Z bosons (Fig. [T]). To achieve signal enhance- 
ment and background reduction, the sum of the trans- 
verse momenta of the two muons and the two leading 
jets, St = PrifJ-i) +PTifJ-2) +PTiji) +PTU2), is also con- 
sidered since the decay products of the SM backgrounds 
are likely to be less energetic than those of the leptoquark 
pair. All combinations of the two muons with the two 
leading jets are taken to calculate four muon-jet invari- 
ant masses M(/i, jet) which are related to the leptoquark 
mass for the signal processes. Together with M(/u, /z)™'" 
and St, these four muon-jet invariant masses M(//,jet) 
define the set of discriminating variables that will be used 
for limit determinations. 

In the case of the iiJ^tH selection, the rate of multi- 
jet events is reduced by requiring the transverse mass 
MTin,^T) = yj2pT{ti)^T [1 - cos A(/)(/i, ^t)] recon- 
structed from the muon and the missing transverse en- 
ergy to be larger than 50 GeV and also by imposing 
]^T > 30 GeV. In order to remove events with mismea- 
sured muon pT, which could lead to large ^t, the az- 
imuthal angle between the missing transverse energy and 
the muon is constrained to be smaller than 3.0 radians. 
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M(n,nr'" (GeV) S^{Gey) 

FIG. 1: [color online] (a) fJ^fJ,jj selection, minimum of the initial and corrected dimuon invariant mass, (b) fi^rjj selection, 
sum of and the transverse momenta of the final state muon and the two jets. Both plots show data (red dots), signal 
(black squares), and background (colored histograms) events. Branching fractions /3 = 1 and (3 = 0.5 and leptoquark masses of 
280 GeV and 200 GeV are assumed for the signal events in (a) and (b), respectively. 



After applying the preceding cuts, the contribution from 
W and Z/j* boson production is simuhaneously nor- 
malized to data with a common scale factor in the re- 
gion 50 GeV < Mrip^, ^t) < HO GeV, which is domi- 
nated by W boson production. At this stage, 5693 data 
events and 5748 ± 395 expected background events re- 
main. When assuming M^q = 200 GeV, the correspond- 
ing signal efficiencies for ^qvq and ^iq^iq events are 22.6% 
and 10.2%, respectively. 

Similar to the dimuon sample, six discriminating 
variables are chosen for the single muon selection. 
Mt{^It ^t) provides a good separation between signal 
and W boson events. The large mass of the hypo- 
thetical leptoquark motivates the choice of the follow- 
ing kinematic variables: the sum of and the trans- 
verse momenta of the final state muon and the two jets, 
St = PT{fJ.)+PT{ji)+PT{j2) + ^T (Fig. [J), the two trans- 
verse masses Mt{^Tt}G^) constructed from the missing 
transverse energy and each of the two highest px jets, 
and the two invariant masses M(/z,jet) derived from the 
muon and the two leading jets. 

For each of the two selections, the six discriminat- 
ing kinematic variables are combined into a multivari- 
ate classifier. A good separation between the leptoquark 
signal and the SM background is obtained with the k- 
Nearest- Neighbors algorithm {kNN), as implemented in 
the TMVA [16| library. The classification relies on the 
comparison of a test event to reference events taken from 
training data sets. The implemented algorithm can be 
interpreted as a generalization of the maximum likeli- 
hood classifier to n dimensions, where n is the number 
of variables used for the discrimination of signal against 
background. During the training phase, the classifica- 
tion of an event as being either signal or background is 



achieved by estimating the local signal-like probability 
density. This is defined as the ratio of the signal events 
over the background plus signal events in the vicinity 
of the tested event, such that the denominator, i.e. the 
number of neighbor events, is equal to the input param- 
eter k. For our analysis, the optimal value for k is found 
to be 50, which results in nearly maximal signal over 
background ratios for reasonable computing times. The 
output of the discriminant takes values between and 1, 
referring to the most background-like events and 1 to 
the most signal-like. 

The training phase is performed for each assumed lep- 
toquark mass in the range between 140 and 320 GeV sep- 
arately (in steps of 20 GeV), and is based on simulated 
signal and background samples for which additional selec- 
tions arc imposed to remove regions with negligible signal 
contribution. For both channels, the corresponding St 
variable is required to be greater than 200 GeV, while 
M{fj,, fj,) and MxifJ., ^t) are required to exceed 100 GeV 
and 110 GeV for the dimuon and the single muon anal- 
yses, respectively. In the case of the fi^xjj selection, 
the shape of the signal event distributions does not vary 
as a function of /3, as the topologies arising from the 
LQLQ — > /i(?/i<7 (with one muon not being reconstructed) 
and fJ-qh'q decay modes are similar. Only one value of /? is 
therefore needed to complete the training for this chan- 
nel. In the case of the dimuon analysis, only the fJ.qfJ-q 
signal events contribute and are considered. The per- 
formance of the training phase for the different assumed 
leptoquark masses is thus also independent of /3. 

In order to avoid overtraining on statistical fluctua- 
tions, the training phase is performed on half of each 
signal and background samples. When using the signal 
and background samples to train the fcAW classifier, alter- 
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(b) D0, L = 1 fb"' 
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FIG. 2: [color online] kNN output variable in the (a) /i/ijj and (b) fi^xjj selections. Both plots show data (red dots), signal 
(filled histograms), and background (full lines) events. Branching fractions /? = 1 and (3 = 0.5 and leptoquark masses of 280 GeV 
and 200 GeV are assumed for the signal events in (a) and (b), respectively. The vertical lines correspond to the borders of the 
bins used in the calculation of the cross section limits. 



natively one event is kept for the training and the next 
one is included in a test sample. Good agreement in 
the kNN distributions between the training and the test 
samples is observed, which demonstrates the absence of 
overtraining. 

Instead of cutting on the kNN variable, we choose a 
strategy to optimize the sensitivity that uses the shape 
of the full kNN distribution. The kNN distributions are 
divided into six bins of variable size which decrease with 
increasing signal efhciency: 0-0.5, 0.5-0.7, 0.7-0.8, 0.8- 
0.9, 0.9-0.95, and 0.95-1. The same binning is used for 
each assumed leptoquark mass and f3 value. The distri- 
butions of the kNN output variable are shown for both 
selections in Fig. [51 The content of each bin is given 
in Tables U and [TT] for the fifijj and ^^Tjj selections, 
respectively. 

The systematic uncertainties on the predicted number 
of background events and on the signal efficiencies are 
studied by varying the efhciencies and resolutions for the 
reconstructed objects and the modeling of both back- 
ground and signal within their uncertainty range [§] . For 
each uncertainty source and for each assumed value for 
P and Mlq, the fcAW output variable is recalculated and 
the deviation observed in each of the kNN bins used for 
the classification is taken as the systematic uncertainty. 
Below we quote the uncertainties obtained when assum- 
ing Mlq — 200 GeV and (3 — 0.5 for the i^J^rij analysis 
and Mlq = 280 GeV and (3 — 1 iov the /i/ijj selection. 
For other assumptions on M^q and (3, the uncertainties 
are similar. 

The dominant systematic uncertainties are listed in Ta- 
ble [TTll For the background they include uncertainties on 



the muon pT resolution, the jet energy scale, the tt cross 
section, and the modeling of jet radiation in W/Z^ jets 
events. The latter is studied by comparing the pT distri- 
bution of the second highest pr jet observed in data with 
the predictions of the ALPGEN and pythia event gener- 
ators when selecting W boson or Z boson events. While 
ALPGEN correctly reproduces the shape observed in data, 
PYTHIA is found to underestimate the rate at large pt- A 
mixture of ALPGEN and pythia events, with a 30% (43%) 
contribution of the latter, gives an acceptable description 
for W+ jets (Z-|-jets) events with a Kolmogorov-Smirnov 
probability corresponding to a la variation. 

Additional uncertainties arise from the efficiencies of 
the muon trigger and the muon identification (2%), of 
the jet reconstruction (0.5%), and from the normaliza- 
tion of the W/Z+ jets background (2%). The rela- 
tive uncertainty on the integrated luminosity is equal to 
6.1% [l3]- The uncertainty on the multijet background 
in the fi^xjj analysis due to the extrapolation into the 
signal region and due to its normalization is estimated to 
be 20%. 

For the signal efficiency and acceptance, the follow- 
ing systematic uncertainties are studied in addition to 
the uncertainties arising from the muon and jet measure- 
ments. The uncertainty due to the PDF choice is deter- 
mined to be 1.6%, using the twenty-eigenvector basis of 
the CTEQ6.1M PDF set The effects of initial and 
final state radiation (ISR and FSR), which might lead 
to the generation of additional jets, are studied by vary- 
ing the PYTHIA parameters controlling the QCD scales 
and the maximal allowed virtualities used in the simula- 
tion of the space-like and time-like parton showers. The 
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TABLE I: Content of each bin of the kN N Yaxiahle after preselection for the /i/xjj analysis. The signal efficiencies for leptoquark 
decay into the fiq^iq state and the expected number of signal events are given, as well as the numbers of events in the data and 
predicted background. The multijet background is negligible. The assumed leptoquark mass is 280 GeV and f3 is taken equal 
to 1. The first uncertainties are statistical and the second systematic. The systematic uncertainty on the integrated luminosity 
is not included. 





n <r kNN <r ^ 


5 kNN <r n 7 


07^ kNN <r08 


Wi^ eiy) + iets 


57 ± 12 ± 07 


0045 ± 0045 ± 0008 




Z/'y*{-^ £+£~) + jets 


48 ± 1 ± 6 


1.00 ± 0.09 ± 0.18 


0.42 ± 0.04 ± 0.03 


tt 


4.4 ± 0.1 ± 0.8 


0.22 ± 0.03 ± 0.08 


0.058 ± 0.018 ± 0.02 


Total background 


53 ± 1 ± 6 


1.2 ± 0.1 ± 0.2 


0.48 ± 0.05 ± 0.04 


Data 


56 





2 




2.94 


2.47 


2.83 


^ signal 


0.69 ± 0.02 ± 0.05 


0.59 ± 0.02 ± 0.02 


0.67 ± 0.02 ± 0.03 


Samples 


0.8 < kNN < 0.9 


0.9 < kNN < 0.95 


0.95 < kNN < 1 


W{~^ lu) +jets 


0.0024 ± 0.0024 ± 0.0005 






zh*{-^ e^e-) +jets 


0.43 ± 0.03 ± 0.15 


0.22 ± 0.03 ± 0.02 


0.31 ± 0.03 ± 0.04 


tt 


0.058 ± 0.016 ± 0.014 


0.031 ± 0.010 ± 0.009 


0.035 ± 0.011 ± 0.007 


Total background 


0.49 ± 0.04 ± 0.16 


0.25 ± 0.03 ± 0.03 


0.34 ± 0.03 ± 0.05 


Data 





1 





<^signal l/cJ 


4.81 


3.97 


18.0 


signal 


1.13 ± 0.02 ± 0.04 


0.94 ± 0.02 ± 0.03 


4.23 ± 0.05 ± 0.22 



TABLE IL Content of each bin of the A;A''A'^ variable after preselection for the fi^xjj analysis. The signal efficiencies for 
leptoquark decay into the fiqi^q and nqfiq states and the expected number of signal events are given, as well as the numbers 
of events in the data and predicted background. The assumed leptoquark mass is 200 GeV and (3 is taken equal to 0.5. The 
first uncertainties are statistical and the second are systematic. The systematic uncertainty on the integrated luminosity is not 
included. 



Samples 


< kNN < 0.5 


0.5 < kNN < 0.7 


0.7 < kNN < 0.8 


W{-^ Iv) +jets 


348 ± 5 ± 44 


6.1 ± 0.5 ± 1.3 


1.8 ± 0.3 ± 0.4 


Z/7*(^ l+r) +jets 


53 ± 1 ± 6 


1.13 ± 0.22 ± 0.15 


0.40 ± 0.07 ± 0.08 


tt 


34.9 ± 0.4 ± 6.5 


2.5 ± 0.1 ± 0.5 


1.0 ± 0.1 ± 0.2 


Multijet 


8.3 ± 0.3 ± 1.7 


0.030 ± 0.017 ± 0.006 


0.020 ± 0.014 ± 0.004 


Total background 


444 ± 5 ± 53 


9.8 ± 0.6 ± 1.7 


3.3 ± 0.3 ± 0.6 


Data 


405 


10 


3 


•^signal 


5.19 


2.35 


1.73 


<Zl^ (%) 


3.57 


1.08 


0.517 


^ ' signal 


8.9 ± 0.1 ± 0.4 


3.9 ± 0.1 ± 0.1 


2.8 ± 0.1 ± 0.2 


Samples 


0.8 < kNN < 0.9 


0.9 < kNN < 0.95 


0.95 < kNN < 1 


Iv) +jets 


1.7 ± 0.3 ± 0.5 


0.31 ± 0.09 ± 0.24 


0.39 ± 0.13 ± 0.18 


zh*{-^ i+e~) +jets 


0.23 ± 0.04 ± 0.01 


0.11 ± 0.02 ± 0.01 


0.064 ± 0.006 ± 0.022 


tt 


0.84 ± 0.05 ± 0.19 


0.26 ± 0.03 ± 0.05 


0.21 ± 0.03 ± 0.06 


Multijet 


0.030 ± 0.017 ± 0.006 


0.020 ± 0.014 ± 0.004 


0.0099 ± 0.0099 ± 0.0020 


Total background 


2.8 ± 0.3 ± 0.5 


0.70 ± 0.10 ± 0.25 


0.67 ± 0.14 ± 0.22 


Data 


4 








<^signal (/o) 


2.28 


1.64 


3.15 


<^signal I/CJ 


0.718 


0.444 


0.600 


signal 


3.7 ± 0.1 ± 0.1 


2.6 ± 0.1 ± 0.1 


5.0 ± 0.1 ± 0.3 



corresponding uncertainty on the signal efficiencies is de- 
termined to be 1.5%. 

For both selections no excess of data over the predicted 
background is observed. Upper hmits at the 95% C.L. 
on the leptoquark production cross section a are cal- 
culated for both selections separately and their combi- 
nation using the method proposed in Ref. [l8| . Since 
both decay modes iJ-qiyq and fJ^qfJ^q contribute in the sin- 



gle muon selection, limits on the product of a and the 
branching fraction cannot be derived and thus the cross 
section limits need to be evaluated for each value of (3 
separately. The six bins (or twelve bins in case of the 
combination) in the kNN discriminant are treated as in- 
dividual channels and their likelihoods are combined with 
correlations of systematic uncertainties taken into ac- 
count. The limits are calculated using the confidence 
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TABLE III: The dominant systematic uncertainties (in %) on the expected number of background events (split between W jZ 
and ti production) and on the leptoquark (LQ) signal efficiency and acceptance for the n^rjj and selections. The 

uncertainty range found for the six bins of each distribution of the kNN discriminant is quoted assuming Mlq = 200 GeV 
and f5 = 0.5 for the fi^rjj selection and Mlq ~ 280 GeV and /3 = 1 for the fijJ.jj analysis. The relative uncertainty on the 
integrated luminosity is 6.1%. 



Uncertainty 


W/Z 


)J.^Tjj channel 
tt 


LQ 


W/Z 


/i^jj channel 
tt 


LQ 


Muon pt resolution 


6-27 


4-23 


0.8-3.2 


2-19 


0-25 


0.1-4.5 


Jet energy scale 


1 - 20 


0-11 


0.1-3.6 


0-11 


0-24 


0.3-4.7 


W/Z+]ets model 


2-77 






4-29 






ti cross section 




18 






18 




ISR/FSR 






1.5 






1.5 


PDF 






1.6 






1.6 



TABLE IV: NLO cross section [T^ for scalar leptoquark pair production in pp collisions at y's = 1.96 TeV using the CTEQ6.1M 
PDF set [l3], observed and expected 95% C.L. cross section limits obtained for the ii^Tjj selection assuming [3 = 0.5 and the 
fj,fj.jj selection assuming P — 1, and observed limits for the combination assuming P = 0.5 or /3 = 1. 







Uprjj channel 


mj 


channel 


Combination 


LQ mass 


(TNLO 


P^ 


0.5 


P 


= 1 


P = 0.5 


P^l 


(GeV) 


(pb) 


CTobs (pb) 


o-cxp (pb) 


o-obs (pb) 


f^oxp (pb) 


o-Qbs (pb) 


Oohs (pb) 


140 


2.38 


0.291 


0.395 


0.0446 


0.0503 


0.127 


0.0426 


160 


1.08 


0.219 


0.204 


0.0392 


0.0357 


0.129 


0.0338 


180 


0.525 


0.144 


0.144 


0.0278 


0.0257 


0.0832 


0.0223 


200 


0.268 


0.0770 


0.106 


0.0262 


0.0206 


0.0551 


0.0220 


220 


0.141 


0.0619 


0.0841 


0.0215 


0.0177 


0.0484 


0.0186 


240 


0.0762 


0.0540 


0.0623 


0.0202 


0.0152 


0.0374 


0.0157 


260 


0.0419 


0.0516 


0.0572 


0.0171 


0.0135 


0.0348 


0.0142 


280 


0.0233 


0.0440 


0.0514 


0.0135 


0.0120 


0.0281 


0.00946 


300 


0.0131 


0.0423 


0.0470 


0.0114 


0.0107 


0.0255 


0.00931 


320 


0.00739 


0.0398 


0.0390 


0.0100 


0.0101 


0.0227 


0.00822 



level CLs — CLs+b/CLb, where CLs+b and CLb are 
the confidence levels for the signal plus b ackg round and 
background only hypotheses, respectively [l8|. 

The observed cross section limits and the expected 
limits are shown in Fig. [3] and Table IIVI together with 
the theoretical prediction for scalar leptoquark pair pro- 
duction calculated at next-to- leading order (NLO) in the 
strong coupling constant using the CTEQ6.1M PDF 
set J^] and the renormalization and factorization scale 
fiR^F = Mlq. Limits for both selections and their com- 
bination are given assuming f3 = 0.5 and /3 = 1. The 
uncertainty band for the cross section prediction shown 
in Fig. [3] reflects the PDF uncertainty IJ] and the vari- 
ation of the factorization and renormalization scale be- 
tween Mlq/^ and 2Mlq, added in quadrature. 

Limits on the leptoquark mass are extracted from the 
intersection of the observed upper bound on the cross 
section with the NLO prediction and also the lower edge 
of its uncertainty band. Combining the fi^xjj and 
fJ-1-J.jj selections and using the central theoretical pre- 
diction, lower bounds on the mass of second genera- 
tion leptoquarks are determined at the 95% C.L. to be 
Mlq > 316 GeV, Mlq > 270 GeV, and Mlq > 185 GeV 
for (3 = 1, P = 0.5, and /3 = 0.1, respectively. Mass 



limits based on the lower edge of the cross section pre- 
diction as well as the expected bounds are listed in Ta- 
ble |Vl Figure m shows the excluded region in the (3 versus 
Mlq parameter space together with the exclusion limits 
obtained for the fj,^Tjj and fJ-^jj selections separately. 
The bound at /3 = 0, where this analysis has no sen- 
sitivity, is given by the DO search in the acoplanar jet 
topology [4| . It is interesting to notice the improvement 
due to the inclusion of the nqfJ,q decay mode in the single 
muon analysis. Therefore, this selection has its maximum 
sensitivity around (3 = 0.6 instead of /3 = 0.5 and a siz- 
able contribution at /3 = 1, where the branching fraction 
for the fJ-qi'q decay mode vanishes. 

In summary, a search for second generation scalar lep- 
toquarks produced in pp collisions at y/s = 1.96 TeV has 
been performed using an integrated luminosity of 1 fb^^. 
Two selections based on the fi^rjj and fJ-fJ^jj final states 
have been carried out. The leptoquark signal was dis- 
criminated from background using a multivariate tech- 
nique based on the kNN algorithm. Both analyses were 
combined to obtain lower limits on the scalar leptoquark 
mass as a function of /3 which exceed 300 GeV at large f3. 
These results improve on previous leptoquark searches 
at the Tevatron [2, They exceed the corresponding 
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TABLE V: Observed and expected mass limits for second generation scalar leptoquarks assuming l3 — 1, (3 = 0.5, and (3 = 0.1. 
The limits are derived using the NLO prediction for the cross section or the lower edge of its uncertainty band. 







Central theory 


Lower edj 


le theory 


13 


M£q (GeV) 


M^^ (GeV) 


Mt4 (GeV) 


ACc? (GeV) 


0.1 


185 


181 


174 


175 


0.5 


270 


272 


259 


263 


1.0 


316 


316 


306 


308 



o. 



10" 



10' 



D0, L = 1 fb"* 

aNLO (Q = MJ 

«jSppe?Lit Observed 
^ ^Sppe?Lit expected 



P=0.5 




P=1 



150 



200 



250 



300 

Mlq (GeV) 



FIG. 3: [color online] Observed and expected 95% C.L. upper 
cross section limits for second generation scalar leptoquark 
pair production assuming j3 = 0.5 or /3 = 1. The NLO pre- 
diction is shown with an uncertainty band reflecting the PDF 
and scale uncertainty. 




nil jj + jj 
(lower edge th.) 

Uli jj + jj 
(central th.) 

Uli jj 

Pr+jets 
(lower edge th.) 

Pr+jets 
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150 
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FIG. 4: The observed 95% C.L. exclusion regions in the Mlq 
versus /3 plane obtained for the fi^rjj and i-L^J.jj selections 
and their combination. The exclusion regions for the two 
separate selections and the solid line for the combination are 
obtained with the cross section prediction reduced by its un- 
certainty. The dashed line is obtained using the nominal NLO 
prediction. The exclusion at vanishing /3 is based on the up- 
dated DO search in the acoplanar jet topology [j. 



previous bounds by 55 GeV at both 13—1 and /? = 0.5, 
and give the most constraining direct limits on second 
generation leptoquarks to date. 
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